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Abstract. - Molecular motors interacting with cytoskeletal filaments undergo peculiar ran- 
dom walks consisting of alternating sequences of directed movements along the filaments and 
diffusive motion in the surrounding solution. An ensemble of motors is studied which interacts 
with a single filament in two and three dimensions. The time evolution of the probability dis- 
tribution for the bound and unbound motors is determined analytically. The diffusion of the 
motors is strongly enhanced parallel to the filament. The analytical expressions are in excellent 
agreement with the results of Monte Carlo simulations. 



Introduction. - Small particles, dispersed in a liquid, constantly undergo random move- 
ments. This was first observed by the botanist Brown and later explained by Einstein and 
Smoluchowski in terms of thermally-excited collisions. The corresponding diffusion coefficient 
depends on the temperature T and on the viscosity of the liquid which reflects the balance 
between thermal fluctuations and viscous dissipation. Such diffusive processes are also ubiq- 
uitous in biological systems M and, in particular, within biological cells 

When the diffusing particles are attracted towards a 2-dimensional surface or a 1-dimen- 
sional filament, they can attain a bound state in which their motion is restricted to this 
low-dimensional manifold. Even though the corresponding diffusion coefficient is expected to 
be small, the reduced dimensionality of the diffusion process can enhance the probability for 
two diffusing particles to collide and, thus, the rate of diffusion-controlled reactions ||. 

Molecular motors can bind to filaments such as microtubules or DNA and then undergo 
1-dimensional diffusion or 'sliding' in the absence of adenosine triphosphate (ATP). This has 
been observed both for cytoskeletal motors at microtubules Q and for RNA polymerase at 
DNA ||. In addition, these motors are able to hydrolyse ATP and to use the released free 
energy in order to change their own conformation. |7) In this way, detailed balance is broken 
and the molecular motors can perform directed walks along the filaments as studied in the 
framework of ratchet models, see, e.g., [|| and the extensive review in 

Since the binding energy, AE, between the motor particle and the filament is necessarily 
finite, the motor particle will again unbind from the filament after a certain walking time, 
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At;,, which is of the order of exp[AE /ksT]. On time scales which are large compared to At;,, 
the motor particle undergoes peculiar random walks which consist of alternating sequences 
of directed movement and unbiased diffusion [H]. Such walks can be studied in the 
context of lattice models which we have recently introduced and studied for open and closed 
compartments using scaling arguments and Monte Carlo (MC) simulations [ jljj . 

In the present work, we consider unbounded geometries, i.e., systems without confining 
walls, for which we are able to obtain analytical solutions both in two and in three dimensions. 
Furthermore, we show that the diffusion of both the bound and the unbound motors is strongly 
enhanced parallel to the filament. For large times, the corresponding diffusion coefficients are 
found to attain anomalously large values in two dimensions and to exhibit large logarithmic 
correction terms in three dimensions. All analytical results are compared with and confirmed 
by extensive MC simulations as shown in Fig. ^ and ^ below. 

2- dimensional case (d = 2). - Consider a discrete time random walk on a square lattice 
with lattice sites labeled by integer coordinates (n, m). At each step, a motor 'particle' has 
probability 1/4 to jump to any of its four nearest neighbor sites. On the line with m = 
which corresponds to the filament, the motion is different: The particle jumps forward from 
(n, 0) to (n + 1, 0) with probability 1 — 7 — kS — \e and backward from (n, 0) to (n — 1, 0) 
with probability i 5, where 7 is the probability to make no jump at all. The latter parameter 
is included in order to incorporate the large difference between the diffusion coefficient in the 
solution and on the filament |]IT| . In addition, the particle leaves the filament and, thus, jumps 
from (n, 0) to (n, ±1) with probability \e where the parameter e is taken to be small. 

As long as the motor is bound to the filament with m = 0, it has the average velocity 
Vi = 1 — 7 — <5 — i e. It also has the probability ^ e to unbind per unit time. Thus, the 
probability that the motor is still bound after t time steps is (1 — | e)* « exp(— i e t). 

As the initial distribution, we take an ensemble of non-interacting particles at n = m = 0. 
The master equation for this dynamics can be solved using Fourier-Laplace (FL) transforms, 
see, e.g., Jl2). The FL transforms for the full probability P and for the probability Pb to be 
bound to the filament are given by 

P(q,r, S )=jr £ ^lyl. PnMt) and P b (r, s) = f] £ -^— P nfi (t) (1) 

and can be related to each other due to the translational invariance of the system parallel to 
the filament. By integration over q, one then obtains the explicit solution 

Pfc ( r ' s ) = —Ti u-i r~~T~ '1 ztt— — — ( 2 ) 

s + [1 — 7J(1 — cosr) + 5 e (cosr — e — iv b smr 

where cosh fi = 2 + 2s — cos r. This result is easily checked for e = (random walk in one 
dimension) and for 7 = 0, S = ^ , and e — 1 (non-biased random walk in two dimensions). 
Likewise, one obtains a closed expression for P(q, r, s) and for the FL transformed distribution 
of the unbound motors which we define via P u b{q, s) = P{q, r, s) — Pb{r, s). 

Fraction of bound motors. - The fraction No(t) = P n ,o(t) of motors bound to the 
filament can be calculated in closed form which leads to 



v„„ , = t £(, + .,.*(, . m = /■'■" 4 » r 4 f^i ,3, 

2ttz J n^/y l + (l-2e)y J Q ir^/y 1 + y 
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Fig. 1 - (a) Fraction No and (b) Average position rib of motors bound to the filament as a function 
of time t. In both figures, the circles and diamonds are MC data for d — 2 and d — 3, respectively, 
and the curves represent the analytical results, see text. For large t, the bound fraction No ~ 1/vt 
in d = 2 (circles) and No ~ 1/t in d = 3 (diamonds). The parameters, which determine the jump 
probabilities of the underlying random walks, have the values 7 = 0, S = 0.6, and e = 0.05. 



where the asymptotic equality holds for small e and large t. For small t, the survival fraction 
N (t) can be expanded in powers of e\/i which leads to Na(t) ~ 1 — 2 ^2 + e 2 t — ie 3 ^ ■ 
Thus, although the motors detach at times ~ they stay close to the filament until 

times t ~ 1/e 2 - For t ^> 1/e 2 , on the other hand, the survival fraction decays to zero as 
No(t) w (1 — + ie d t i )/Vvr e 2 t. The complete time evolution of No(t), as obtained by 
numerical evaluation of the exact integral in (3), is displayed in Fig. [l](a). We have also 
determined this quantity in a completely different way using Monte Carlo (MC) simulations. 
Inspection of Fig. |](a) shows that the results of these two approaches agree very well. 

Average position on the filament. - The first moment of the probability distribution P n fi 
for the motors on the filament is given by 

Wl(t)sE „p„, o(()= ^r^iz^ (4) 

„ Tre Jo VV I 1 + V) 

as follows from an expansion of (2) to first order in r. For small t, this leads to N%(t) w 

v b t [1 - (8/3-y/n)eVi H ]. In the same limit, the average position = N±/No and the 

average velocity tJf, = dno/dt behave as rib ~ Vbt(l — ^ey/t/^/n) and Ub w Vb(l — ^^); 
respectively. For large t, the integral in (Q) leads to Nx(t) w ^=^) and thus to 



n b w (1 =) and v b « -^oW v b w — ^ • (5) 



The complete time evolution of the average position n& is displayed in Fig. |l](b). Again, we 
compare rib as obtained from the integrals in (3) and (4) with MC data and find very good 
agreement. Note that the average position Ub of the bound motors behaves as ~ \ft for large 
t. As shown below, the same time dependence applies to the average position of all motors as 
previously predicted by scaling arguments (lQ, O 
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Variance and diffusion coefficient on the filament. - The second moment of the proba- 
bility distribution P Ui o for the bound motors is 

^— ' 7re 4 7 (1 + x^ar/^ 

as obtained from an expansion of (0) to second order in r and for small e . For small t, this leads 
toJV 2 (t) » (l- 7 )i+v b 2 i 2 (l-^ ). The positional variance, defined via An 2 = ^§y-^l§y, 
i.e., with respect to the conditional probability P n ,o(t)/No(t), then behaves as 

2m 2 et 5 / 2 

An|« (l-7)t+ — for small* . (7) 

15 V 71 " 

The relative variance An 2 /n 2 « ^ ^= is small since t -C 1/e 2 - We now define a iime- 
dependent diffusion coefficient via Db(t) = \ dAn 2 /di which exhibits the small-t behavior 

A(t)«i(l-7) + ^«6 2 e^ 2 • (8) 

In addition, the behavior for large t is found to be 

Anl* V -l(4-n-^)t and D b {t) » |J(4 - - ^) . (9) 

Note that the diffusion constant Db{oo) ~ v 2 /2e 2 is quite large compared to the value Df,(0) = 
2 (1 — 7) as appropriate for the 1-dimensional diffusion of a bound motor. In fact, one has 
L>b(t) ~ f 2 /e 2 as soon as t <; 1/e 2 , see (||). This broadening occurs since the unbound motors 
lag behind the bound ones. The time-dependence of Art 2 is shown in Fig. 0. 

Probability distribution on the filament. At large times, N2/N0 ~ (Ni/Nq) 2 which 
indicates that the probability distribution attains a scaling form. It follows from (^) that 
Pb(r, s) w l/(s — iwfcr + e-^/s) for small r and s. The inverse Fourier transform gives P n .o{s) = 
for n < 0. For n > 0, we close the contour along the negative real axis. After the substitution 
s = —u 2 , the u- integral can be taken over the whole real axis, at the expense of a factor | . 
The integral is then Gaussian and leads to 



PnfiW ~ o ; / 7 TTTJ ex P 

2^/irvb(Vbt — n)*' z 



4:V b (v b t 



for < n < v b t , (10) 



which vanishes exponentially in l/(u{,£ — n) as n approaches v&i from below which shows that, 
for large i, all motors unbind before they reach the position n = v b t. 

Properties of the unbound motors. Eventually all motors will become unbound and 
diffuse in the surrounding solution. We now discuss how this behavior is affected by the 
filament. The FL transformed distribution, P u b(q,r, s), of the unbound motors behaves as 

P ub (q, r, s) « : P^-7=r, 2 , a \ ( n ) 

for small q, r, s with q ~ y/s. Using this asymptotic expression , we can calculate the moments 
Nl ib (t) = EnEm/o 11 ' 1 ^'-"' f° r tne distribution of the unbound motors. The fraction of 
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Fig. 2 - Positional variance An| of bound motors as a function of time t. The MC data for d — 2 
(cirles) and d = 3 (diamonds) are again compared with the analytical results (full and dashed curve). 
For large t, the diffusion coefficient Db(t) — ~ dAnl/dt attains the anomalously large value Db(oo) = 

^t-(4 — 7r) in d = 2 (circles) whereas Dt(t) ~ | + ^7-5 ^ in d = 3 (diamonds) with the usual 
diffusion constant £)j,(oo) = A but a large correction term ~ \nt/t. The jump probabilities have the 
same values as in Fig. M. 



unbound motors, N^ b , satisfies Ntf b (t) = 1 — No(t). The first moment iV" b determines the 
average displacement and velocity, n u b and ivt,, of the unbound motors parallel to the filament 
via n u b = Ni b /NQ b and v u b = dn u b/dt. For small t, n u b f=s §f&i(l — ^ ey/ni) and v U 6 w 
§i>b(l — § e\fm), while for large t, 

H*«*^[l-£-l)-L=] and ^ = _!* [1 + O(i)]. (12) 



Although each individual unbound motor has, on average, zero velocity, the velocity v u b is 
non-zero, since the cloud of unbound motors moves via repeated interactions with the filament. 
The average velocity of both bound and unbound motors is given by NoVb + NQ b v u b- For large t, 
this latter quantity is dominated by the unbound term NQ b v u b which becomes asymptotically 
equal to NoVb in agreement with scaling arguments uM- 

For small and large t, the positional variance An^behaves as 

Ar&« « 6 V±(l-i^) and An 2 ub *2(n-2)v b 2 t/ne 2 , (13) 

respectively. The corresponding diffusion constant is given by Dm (00) = (tt — 2)vb 2 /tt e 2 . The 
relation Dm (00) ~ Vb 2 / e 2 implies that the parallel diffusion in the bulk is strongly enhanced 
by the presence of the filament. 

The motor displacement perpendicular to the filament has the variance 

with the asymptotic behavior 

9et 3 / 2 1 1 

Am 2 « — -=■ and Am 2 w - t (1 - 2—==) (15) 
3v7T 2 eyirt 
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for small and large t, respectively. Thus, the perpendicular diffusion is initially suppressed but 
the corresponding diffusion coefficient eventually attains its unperturbed value D±(oo) = |. 
In real space, the probability distribution as given by (O) becomes 



e(en + 2\m\v b ) 
= o 7=77. I _ v]3 /2 ex P 



(en + 2|m|ufe) 



4:Vb(vbt - n) 



for 0<n<v b t (m^O). (16) 



3-dimensional case (d — 3). - Now consider the same kind of random walk but on a cubic 
lattice. The jump probabilities per time step are g away from the filament, 1 — 7 — \ 8 — |e 
and \ 5 in the forward and backward direction, respectively, on the filament, and g e for each 
of the four unbinding directions, where 7 is again the probability to make no step at all. 
The average velocity within the bound state is Vb = 1 — 7 — S — | e , and the probability 
to stay in such a bound state for t successive steps is exp[— 2e t/3]. We now denote the 
perpendicular coordinate by m = (mi, 7712) and the perpendicular momentum by q = (gi, q 2 ). 
After integration over qi and (72, we obtain 

Ph (rs) = 3/ ( r ' s ) (17) 

e+[3(l-e)s+i(e-3J)(e-^-l)-i(6-67-3^-5e)(e»--l)]/(r,s) V ; 

with 

77 \- dgi dg 2 1 ^ , 1SJ , 

i(r, s) = / = K(m), (18) 

Jq 2tt 2tt 3 + 3s — cos r — cos qi — cos g 2 tt 

where K(m) = J* dtp/ \fl — m sin 2 <f> is a complete elliptic integral at m = 4/(3 + 3s — cosr) 2 . 

Using the explicit expressions for the FL transform of the probability distributions, one can 
now compute all quantities of interest for the 3-dimensional case. Since these computations 
are somewhat tedious, we will only present the results and describe the details elsewhere fl3f . 

The statistical properties of the bound motors are displayed in Figs. [I] and || where the 
fraction Nq, the average position fib, and the variance An 2 are shown both for the 2- and 
the 3-dimensional case. For large t, the displacement of the unbound motors behaves as 
n u b ~ 2^e"( mi + 7b)- Again, this is equal to the asymptotic displacement averaged over both 
bound and unbound motors, and the corresponding velocity is given by VbNo(t) ~ 3wb/(27rei) 
as predicted by scaling arguments jyj. 

The parallel variances of the bound and unbound motor positions exhibit logarithmic 
corrections and behave as An 2 (t) w | +*(t) and An 2 f) (t) w | + for large f with ^(t) = 

2^.2^ 2 [(lni + 7b) 2 — ^jj-], tq = 3/16, and Euler's constant — 0.5778. The corresponding 
time-dependent diffusion coefficients behave as Db{t) ~ \ + and D^(t) w i + 

with ^'(t) = 2" 6 2 ; [Ini + 7b]- At the crossover time t ~ 1/e, the logarithmic terms are 
~ ln(l/ e )/ e and, thus, are large compared to the value | which is the diffusion coefficient of 
the motors in the absence of the filament. For displacements in the directions perpendicular 
to the filament, one finds the variances Am 2 (£) = Am|(t) ~ | — 577 + 7b) f° r large t. 

Filament with several proto filaments. Microtubules consist of 13 parallel protofila- 
ments, on each of which motors can run. To incorporate this in our model, we assume that 
the pinning line at m = has k internal states, and k may be equal to 13. The average 
occupation of each of these states is denoted by , j = 1, .., k. There is a small probability, 
\ £, that a motor goes from protofilament j to j + 1 in one time step, and similarly for going 
to protofilament j — 1. To take into account the cylindrical structure of the microtubule, 
we identify j = k with j = 0. We assume that, after unbinding, the motor may enter the 
surrounding solution at a randomly chosen neighboring site of our lattice. 
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In this model, the total fraction of motors at position n along the tubule is P nt o = 
2Zj—i P-nO- The motion on the individual protofilaments can again be described by a mas- 
ter equation. To solve it, we introduce the Fourier transform P£o — Sj=i PnV'''" with 
u> = 2ir£/k with £ — 0, .., k — 1. At the initial time t = 0, all motors are taken to be on 
the filament position n = 0, m = 0, and on the protofilament j = 0. We then find for the 
asymmetric modes with u =/= 

Pb " (r ' S) = 5 + 1- (1-7-15 - |e)e- - 7 - \ Se~ ir + ((1 - cose;) (19) 
which, in this simple model, has decoupled from P^ =0 (r, s). 



The bound fractions on the individual protofilaments are obtained from (19) with r = 
which leads to 

N j (t) ee $>£ (f) = iiVb(t) + ~ ^ e_yW e-^e-^ 1 - 008 -)* (20) 

The first term describes the symmetric distribution of the motors over the k protofilaments 
and decays algebraically, see Fig. |l](a). The second term represents the asymmetry arising 
from the initial distribution. There are two decay mechanisms of this asymmetry. The term 
exp(— |et) expresses that unbinding and subsequent rebinding to a randomly chosen protofil- 
aments restores the symmetry. The factor exp[— £(1 — cosw)i] expresses that hopping to 
neighboring protofilaments also restores the symmetry. The asymmetry between the average 
occupation of the various protofilaments thus disappears after some transient time. 

In summary, we have studied the walks of molecular motors which arise from many diffu- 
sional encounters with a single filament. In the absence of confining walls or boundaries, we 
have been able to analytically calculate the statistical properties of these walks both for two 
and for three dimensions. The analytical expressions were found to be in excellent agreement 
with the results of MC simulations. The diffusion of both bound and unbound motors was 
shown to be strongly enhanced parallel to the filament. This enhanced diffusion could be 
used, e.g., in order to increase the rate of diffusion-controlled chemical reactions. 
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